We perform reproducibility analyses for ANTsPyMM based on traveling subject data.
SRPBS (T1, rsfMRI): Tanaka Saori, C., et al. “SRPBS Traveling Subject MRI Dataset.” (No Title) (2019).
Tong (T1, DTI): Tong, Qiqi, et al. “Multicenter dataset of multi-shell diffusion MRI in healthy traveling adults with identical settings.” Scientific Data 7.1 (2020): 157.
Hawco (T1, rsfMRI, DTI): Hawco, Colin, et al. “A longitudinal human phantom reliability study of multi-center T1-weighted, DTI, and resting state fMRI data.” Psychiatry Research: Neuroimaging 282 (2018): 134-142.
Each dataset is analyzed by distributed computing with ANTsPyMM as in
src/hawco_processing.
Traveling subject studies involve scanning the same subjects on multiple MRI scanners at different locations. These studies help in assessing consistency and/or agreement of image quantification where the only variables are the machines themselves. This is crucial for understanding power in multi-site studies of natural history or intervention and for ensuring that the observed changes in brain structure or function are due to actual physiological changes rather than variations in the imaging process itself.
In this study, we employ traveling cohort data to assess the agreement of imaging derived phentoypes (IDPs) pooled across multiple sites for the purposes of statistical inference. These data will establish expectations of repeatability for T1w, DTI and rsfMRI as measured by ANTsPyMM processing. Thus, we use these data to characterize the consistency and reliability of these tools when applied to data that has known systematic biases due to site and scanner differences. The results confirm that findings and conclusions drawn from ANTsPyMM are reliable and not overwhelmed by scanner-specific differences or inconsistencies. This knowledge is critical for a foundational framework such as ANTsX upon which scientific studies, machine learning platforms and other methodological comparisons are based.
we investigate three traveling cohort datasets and ANTsPyMM reproducibility therein. these cohorts represent variability in both MRI manufacturer and MRI model (high variability) that would exceed standard (within-scanner, within-site) test-retest analysis. results therefore provide a lower-bound on reliability; i.e. within-site (e.g. longitudinal) studies would be expected to have higher reliability in general.
a traveling cohort : 9 healthy subjects travel to 12 sites to be imaged
of the 12 sites, 9 have consistently available T1w and rsfMRI in 6 subjects
a traveling cohort : 3 healthy subjects travel to 4 sites to be imaged
available at figshare in BIDS format (here).
a traveling cohort : 4 healthy subjects travel to 6 sites to be imaged
“Data are available from four healthy male participants, aged 34 to 59 during study start (aged 38 to 63 at study end). Participants had no history of psychiatric, neurological, or major chronic illness.” Note: true ages are unavailable but were estimated with antspynet.brain_age.
29 subjects/visits have complete data ie all 3 modalities. We select only the baseline data for each site.
these data are also longitudinal with followup 1 to 3 years duration from baseline.
These data enable us to investigate the reliability of our imaging-derived phenotypes (IDPs) computed with ANTsPyMM (latest version)
we use the intraclass correlation coefficient (ICC) to assess consistency or reproducibility of the quantitative IDPs
this analysis focuses on the baseline data ( i.e. we do not investigate the longitudinal effects )
site CMP in the Hawco cohort exists only at year 3 and – as such – does not contribute to baseline data.
We find that (1) the majority of IDPs show superior reliability and (2) joint reduction of IDPs with SiMLR improves reliability further. ANTsPyMM IDPs derived from the same subjects imaged at different sites with MRI from various manufacturers show overall high reliability. This provides empirical evidence that multiple modality MRI may be used to derive quantitative phenotypes on which predictive models may be based.
FIXME - move this statement: Statistical control for site effects should still be applied at the population level.
see this paper 10.1016/j.jcm.2016.02.012 for discussion of ICC
Cicchetti (1994) gives the following often quoted guidelines for interpretation for kappa or ICC inter-rater agreement measures:
A different guideline is given by Koo and Li (2016):
see this analysis of freesurfer on T1w showing values from 0.81 - 0.88
FIXME
T1Hier_resnetGrade is a deep learning based method that
accurately predicts image quality in multi-site data.
Values range from 0 ( unusable ) to 3 ( best ) quality.
##
## ATT ATV CMH COI HKH KPM MRC MRP
## 6 6 4 6 6 6 3 1
## ses-c09r1 ses-c10r1 ses-c10r2 ses-c10r3 SWA UTO YC1 YC2
## 3 3 3 3 6 6 6 6
## ZHH ZHP
## 3 1
## Warning in RColorBrewer::brewer.pal(n, pal): n too large, allowed maximum for palette Set3 is 12
## Returning the palette you asked for with that many colors
## Warning: Removed 70 rows containing missing values (`geom_point()`).
What is the variability of the measurement if we control for year alone?
mdl=(lm( outcome ~ Year + Subject , data=dd ))
visreg::visreg( mdl, "Subject", gg=TRUE ) +
scale_color_brewer(palette = brewpal )+theme(axis.text.x = element_text(angle = 45, vjust = 0.5, hjust=1))
What is the reproducibility of the measurement if we control for age and scanner?
###############################################
mdl=( lm( outcome ~ age_BL + Subject, data=dd ) )
visreg::visreg( mdl, "Subject", gg=TRUE ) +
scale_color_brewer(palette = brewpal )+theme(axis.text.x = element_text(angle = 45, vjust = 0.5, hjust=1))
# grid.arrange( grobs = visreg::visreg( mdl, gg=TRUE ), ncol=1, main='reproducibility' )
What is the reproducibility of the measurement if we control for site?
mdl=(lm( outcome ~ age_BL + Subject + Site , data=dd ))
visreg::visreg( mdl, "Subject", gg=TRUE ) +
scale_color_brewer(palette = brewpal )+theme(axis.text.x = element_text(angle = 45, vjust = 0.5, hjust=1))
# grid.arrange( grobs = visreg::visreg( mdl, gg=TRUE ), ncol=1, main='reproducibility' )
What is the effect of age if we control for the best confounds?
library(lme4)
mdl=(lm( outcome ~ age_BL + T1Hier_resnetGrade +(Subject), data=dd ))
# mdl=(lm( outcome ~ age +(Subject), data=dd ))
visreg::visreg( mdl, 'age_BL', gg=TRUE, main="..." ) +
scale_color_brewer(palette = brewpal )
We elect to use ICC(2,k) (Two-Way Random, Absolute Agreement, Average Measures): This version of ICC is used when the measurement devices are considered random samples from a larger population of raters and we prefer findings to generalize. ICC(2,k) is appropriate when assessing the reliability of subject-specific measurements across sites. The ‘k’ form of ICC means that the reliability is based on the average of multiple measurements (in this case, multiple imaging sites), which generally provides a more robust and stable estimate of reliability.
We study the consistency of MRI IDPs collected from different imaging locations or scanners and processed by ANTsPyMM. These multiple measurements represent a realistic range of conditions under which MRI IDPs are obtained in practice. Such data is typically interpreted via population-level statistical summaries which often control for both site and subject level random effects. In this context, the ICC(2,k) model is the appropriate choice as it considers both between-subject variability and measurement variability across different scanners. In contrast, ICC(1,k) considers scanners as random effects while ICC(3,k) considers subjects as random effects. ICC(2,k) considers both which is analogous to standard population-level statistical modeling where both site and subject variables may be treated as random effects.
print( reli_stats( "outcome", "Site", "Subject", data=dd ) )
##
## Coefficient of Variation (%): -1.68
## Standard Error of Measurement (SEM): 0.0376
## Standard Error of the Estimate (SEE): 0.285
## Standard Error of Prediction (SEP): 0.536
##
## Intraclass Correlation Coefficients with 95 % C.I.
## Model Measures Type ICC Lower CI Upper CI
## 1 one-way random Agreement ICC1 -0.22072 -2.408e-01 -0.03434
## 2 two-way random Agreement ICC2 0.01446 -1.636e-04 0.14249
## 3 two-way fixed Consistency ICC3 0.39352 4.142e-02 0.87839
## 4 one-way random Avg. Agreement ICC1k -9.42215 -3.276e+01 -0.19901
## 5 two-way random Avg. Agreement ICC2k 0.06833 -8.187e-04 0.45380
## 6 two-way fixed Avg. Consistency ICC3k 0.76439 1.777e-01 0.97306
print( reli_aov( "outcome", "Site", "Subject", data=dd ) )
##
## Coefficient of Variation (%): -2.18
## Standard Error of Measurement (SEM): 0.0493
## Standard Error of the Estimate (SEE): 0.125
## Standard Error of Prediction (SEP): 0.181
##
## Intraclass Correlation Coefficients with 95 % C.I.
## Model Measures Type ICC Lower CI Upper CI
## 1 one-way random Agreement ICC1 0.5177 0.1133 0.9608
## 2 two-way random Agreement ICC2 0.6879 0.2747 0.9775
## 3 two-way fixed Consistency ICC3 0.9044 0.6628 0.9947
## 4 one-way random Avg. Agreement ICC1k 0.8430 0.3899 0.9919
## 5 two-way random Avg. Agreement ICC2k 0.9168 0.6544 0.9954
## 6 two-way fixed Avg. Consistency ICC3k 0.9793 0.9076 0.9989
##
## Coefficient of Variation (%): 12.7
## Standard Error of Measurement (SEM): 0.13
## Standard Error of the Estimate (SEE): 0.182
## Standard Error of Prediction (SEP): 0.374
##
## Intraclass Correlation Coefficients with 95 % C.I.
## Model Measures Type ICC Lower CI Upper CI
## 1 one-way random Agreement ICC1 0.2720 -0.0234567 0.8272
## 2 two-way random Agreement ICC2 0.2835 -0.0007801 0.8278
## 3 two-way fixed Consistency ICC3 0.3081 -0.0153657 0.8448
## 4 one-way random Avg. Agreement ICC1k 0.6513 -0.1294271 0.9599
## 5 two-way random Avg. Agreement ICC2k 0.6643 -0.0039129 0.9601
## 6 two-way fixed Avg. Consistency ICC3k 0.6900 -0.0818596 0.9646
## [1] "rsfMRI_fcnxpro122_falffPoint189_STLRAVG"
##
## Coefficient of Variation (%): 13.6
## Standard Error of Measurement (SEM): 0.062
## Standard Error of the Estimate (SEE): 0.205
## Standard Error of Prediction (SEP): 0.297
##
## Intraclass Correlation Coefficients with 95 % C.I.
## Model Measures Type ICC Lower CI Upper CI
## 1 one-way random Agreement ICC1 0.3287 0.01276 0.8528
## 2 two-way random Agreement ICC2 0.4040 0.09556 0.8618
## 3 two-way fixed Consistency ICC3 0.9199 0.75890 0.9903
## 4 one-way random Avg. Agreement ICC1k 0.7100 0.06071 0.9666
## 5 two-way random Avg. Agreement ICC2k 0.7722 0.34566 0.9689
## 6 two-way fixed Avg. Consistency ICC3k 0.9829 0.94026 0.9980
## [1] "rsfMRI_fcnxpro122_alffPoint122_InsLRAVG"
##
## Coefficient of Variation (%): 2.7
## Standard Error of Measurement (SEM): 0.00274
## Standard Error of the Estimate (SEE): 0.0136
## Standard Error of Prediction (SEP): 0.0197
##
## Intraclass Correlation Coefficients with 95 % C.I.
## Model Measures Type ICC Lower CI Upper CI
## 1 one-way random Agreement ICC1 0.2591 0.0845 0.6175
## 2 two-way random Agreement ICC2 0.3133 0.1153 0.6539
## 3 two-way fixed Consistency ICC3 0.9167 0.8256 0.9766
## 4 one-way random Avg. Agreement ICC1k 0.7589 0.4537 0.9356
## 5 two-way random Avg. Agreement ICC2k 0.8042 0.5398 0.9444
## 6 two-way fixed Avg. Consistency ICC3k 0.9900 0.9771 0.9973
## [1] "DTI_mean_fa_LRAVG_pars_orbitalis"
## boundary (singular) fit: see help('isSingular')
##
## Coefficient of Variation (%): 17.2
## Standard Error of Measurement (SEM): 0.145
## Standard Error of the Estimate (SEE): 0
## Standard Error of Prediction (SEP): 0.364
##
## Intraclass Correlation Coefficients with 95 % C.I.
## Model Measures Type ICC Lower CI Upper CI
## 1 one-way random Agreement ICC1 0 -0.1604 0.6061
## 2 two-way random Agreement ICC2 0 -0.1665 0.6077
## 3 two-way fixed Consistency ICC3 0 -0.1665 0.6077
## 4 one-way random Avg. Agreement ICC1k 0 -2.2389 0.8850
## 5 two-way random Avg. Agreement ICC2k 0 -2.4903 0.8856
## 6 two-way fixed Avg. Consistency ICC3k 0 -2.4903 0.8856
## [1] "rsfMRI_fcnxpro122_falffPoint248_StriateLRAVG"
## boundary (singular) fit: see help('isSingular')
##
## Coefficient of Variation (%): 20.4
## Standard Error of Measurement (SEM): 0.185
## Standard Error of the Estimate (SEE): 0.383
## Standard Error of Prediction (SEP): 0.571
##
## Intraclass Correlation Coefficients with 95 % C.I.
## Model Measures Type ICC Lower CI Upper CI
## 1 one-way random Agreement ICC1 0.8177 0.5548 0.9759
## 2 two-way random Agreement ICC2 0.8177 0.5548 0.9759
## 3 two-way fixed Consistency ICC3 0.8177 0.5332 0.9761
## 4 one-way random Avg. Agreement ICC1k 0.9573 0.8617 0.9951
## 5 two-way random Avg. Agreement ICC2k 0.9573 0.8617 0.9951
## 6 two-way fixed Avg. Consistency ICC3k 0.9573 0.8510 0.9951
## [1] "rsfMRI_fcnxpro122_alffPoint025_PFClLRAVG"
##
## Coefficient of Variation (%): 48.5
## Standard Error of Measurement (SEM): 1.01
## Standard Error of the Estimate (SEE): 2.24
## Standard Error of Prediction (SEP): 3.5
##
## Intraclass Correlation Coefficients with 95 % C.I.
## Model Measures Type ICC Lower CI Upper CI
## 1 one-way random Agreement ICC1 0.5047 0.1499 0.9122
## 2 two-way random Agreement ICC2 0.5280 0.1882 0.9138
## 3 two-way fixed Consistency ICC3 0.6903 0.3315 0.9546
## 4 one-way random Avg. Agreement ICC1k 0.8359 0.4686 0.9811
## 5 two-way random Avg. Agreement ICC2k 0.8483 0.5369 0.9815
## 6 two-way fixed Avg. Consistency ICC3k 0.9177 0.7126 0.9906
## [1] "rsfMRI_fcnxpro122_alffPoint018_PFCdLRAVG"
## boundary (singular) fit: see help('isSingular')
##
## Coefficient of Variation (%): 20.4
## Standard Error of Measurement (SEM): 0.287
## Standard Error of the Estimate (SEE): 0
## Standard Error of Prediction (SEP): 1.1
##
## Intraclass Correlation Coefficients with 95 % C.I.
## Model Measures Type ICC Lower CI Upper CI
## 1 one-way random Agreement ICC1 -0.1402 -0.21392 0.3195
## 2 two-way random Agreement ICC2 0.0000 -0.05814 0.3736
## 3 two-way fixed Consistency ICC3 0.0000 -0.16645 0.6077
## 4 one-way random Avg. Agreement ICC1k -1.5970 -7.41144 0.7012
## 5 two-way random Avg. Agreement ICC2k 0.0000 -0.37880 0.7489
## 6 two-way fixed Avg. Consistency ICC3k 0.0000 -2.49029 0.8856
## [1] "rsfMRI_fcnxpro122_perafPoint115_OFCLRAVG"
## boundary (singular) fit: see help('isSingular')
##
## Coefficient of Variation (%): 16.4
## Standard Error of Measurement (SEM): 0.165
## Standard Error of the Estimate (SEE): 0.255
## Standard Error of Prediction (SEP): 0.473
##
## Intraclass Correlation Coefficients with 95 % C.I.
## Model Measures Type ICC Lower CI Upper CI
## 1 one-way random Agreement ICC1 0.4107 0.07140 0.8837
## 2 two-way random Agreement ICC2 0.4107 0.07052 0.8837
## 3 two-way fixed Consistency ICC3 0.4107 0.05388 0.8843
## 4 one-way random Avg. Agreement ICC1k 0.7770 0.27770 0.9743
## 5 two-way random Avg. Agreement ICC2k 0.7770 0.27502 0.9744
## 6 two-way fixed Avg. Consistency ICC3k 0.7770 0.22163 0.9745
## [1] "rsfMRI_fcnxpro122_perafPoint086_PFCdLRAVG"
##
## Coefficient of Variation (%): -3.38
## Standard Error of Measurement (SEM): 0.0915
## Standard Error of the Estimate (SEE): 0.23
## Standard Error of Prediction (SEP): 0.343
##
## Intraclass Correlation Coefficients with 95 % C.I.
## Model Measures Type ICC Lower CI Upper CI
## 1 one-way random Agreement ICC1 0.5220 0.1659 0.9169
## 2 two-way random Agreement ICC2 0.5526 0.1911 0.9196
## 3 two-way fixed Consistency ICC3 0.8126 0.5237 0.9753
## 4 one-way random Avg. Agreement ICC1k 0.8452 0.4986 0.9822
## 5 two-way random Avg. Agreement ICC2k 0.8606 0.5416 0.9828
## 6 two-way fixed Avg. Consistency ICC3k 0.9559 0.8461 0.9950
## [1] "rsfMRI_fcnxpro122_VisPeri_2_SomMotB"
## boundary (singular) fit: see help('isSingular')
##
## Coefficient of Variation (%): 15.7
## Standard Error of Measurement (SEM): 0.163
## Standard Error of the Estimate (SEE): 0
## Standard Error of Prediction (SEP): 0.501
##
## Intraclass Correlation Coefficients with 95 % C.I.
## Model Measures Type ICC Lower CI Upper CI
## 1 one-way random Agreement ICC1 -0.07904 -0.19171 0.4742
## 2 two-way random Agreement ICC2 0.00000 -0.09943 0.4954
## 3 two-way fixed Consistency ICC3 0.00000 -0.16645 0.6077
## 4 one-way random Avg. Agreement ICC1k -0.57792 -4.11069 0.8185
## 5 two-way random Avg. Agreement ICC2k 0.00000 -0.82539 0.8307
## 6 two-way fixed Avg. Consistency ICC3k 0.00000 -2.49029 0.8856
## [1] "rsfMRI_fcnxpro122_perafPoint099_PHCLRAVG"
Site-wise reliability
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in qf(1 - alpha, v, n_id - 1): qbeta(a, *) =: x0 with |pbeta(x0,*) -
## alpha| = 0.049799 is not accurate
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in qf(1 - alpha, v, n_id - 1): qbeta(a, *) =: x0 with |pbeta(x0,*) -
## alpha| = 0.042171 is not accurate
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
## anat icc
## 40 fcnxpro122_alffPoint204_21LRAVG 0.1363677
## 1181 fcnxpro122_perafPoint180_6LRAVG 0.8376975
higher scores are better
a few sites are lower quality
a few subjects exhibit lower quality (consistently)
##
## DTI_fa_gm DTI_fa_jhu rsfMRI_alff rsfMRI_cnx rsfMRI_dfn
## 60 19 223 205 51
## rsfMRI_dfn_att rsfMRI_dfn_tsk rsfMRI_falff rsfMRI_peraf t1_cit
## 24 18 129 193 16
## t1_nbm_mtl t1_thk_ctx t1_vol_cbl t1_vol_ctx
## 13 31 13 31
| modality_st | icc | sd | |
|---|---|---|---|
| 13 | t1_vol_cbl | 0.9767645 | 0.0132082 |
| 11 | t1_nbm_mtl | 0.9728687 | 0.0345903 |
| 10 | t1_cit | 0.9707994 | 0.0194247 |
| 14 | t1_vol_ctx | 0.9681434 | 0.0192710 |
| 12 | t1_thk_ctx | 0.9501737 | 0.0274131 |
| 5 | rsfMRI_dfn | 0.9102535 | 0.0480177 |
| 4 | rsfMRI_cnx | 0.9085471 | 0.0544152 |
| 6 | rsfMRI_dfn_att | 0.8966800 | 0.0588684 |
| 7 | rsfMRI_dfn_tsk | 0.8953134 | 0.0672987 |
| 2 | DTI_fa_jhu | 0.8870767 | 0.0941680 |
| 1 | DTI_fa_gm | 0.8850481 | 0.1256671 |
| 3 | rsfMRI_alff | 0.7024170 | 0.2515456 |
| 9 | rsfMRI_peraf | 0.6766020 | 0.2589057 |
| 8 | rsfMRI_falff | 0.5027130 | 0.2592247 |
train simlr on each of top quality sites
association of T1-derived IDPs and rsfMRI IDPs
default parameters for SiMLR (regression and ICA)
test on each other site
this lets us look at ICC in the simlr space
## [1] "USE SITE: CMH"
## Warning in sqrt(ICC3 * (1 - ICC3)): NaNs produced
rsfMRI measurements are inter and intra network connectivity between the canonical functional networks
default mode
salience
frontoparietal task control
…